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Abstract 

In the present paper we give the general solution of the functional equation 

(f(x) + g(y) + h(z)) 2 = F(x) + G(y) + H(z), x+y+z=0 

which is related to the exact factorized ground-state wave function for the 
quantum one- dimensional problem of three different particles with pairwise 
interaction. 
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Functional equations connecting several functions and admitting a general 
analytic solution have recently attracted the attention of many mathemati- 
cians as well as physicists (for recent results, see, for example, [BC 1990], 
[BP 1993], [BK 1993], [BFV 1994]). 

In modern mathematical physics such equations arise in connection with 
the integrable systems of classical and quantum mechanics (see, for example, 
reviews [Pe 1990], [OP 1983]). 

In the present paper (a previous version of which appeared as [BP 1993]), 
we investigate one such equation. Namely, we investigate the functional 
equation connecting six unknown functions 

(f(x) + g(y) + h{z)f = F{x) + G{y) + H(z), x + y + z = 0, (1) 

which generalizes the well-known Frobenius-Stickelberger equation [FS 1880] 
and is related to the exact factorized ground-state wave function for the 
quantum one-dimensional problem of three different particles with pairwise 
interaction. We give the general non-degenerate solution of this equation. 

1. Let us recall first that an analogous (but simpler) equation for the spe- 
cial case of three identical particles was considered earlier by B. Sutherland 
[Su 1975] and F. Calogero [Ca 1975]. Namely, in the paper [Su 1975] the one- 
dimensional many-body problem of n identical particles with pair interaction 
was considered, whose exact ground-state wave function ^0(^1, ^2, • • • , x n ) is 
factorized 

ty (x 1 ,x 2 ,...,x n ) = Y[ip(x j -x k ). (2) 

j<k 

It was shown that the logarithmic derivative of ^(x) 

f(x) = ^'(x)/^(x) (3) 
must satisfy the functional equation 

/(*) f(v) + f(v) /(*) + /(*) /(*) = F(x) + F(y) + F(z), 

x + y + z = o^ { ) 

where f(x) (F(x), respectively) is an odd (even) function 

/(-*) = -/(*), F(-x) = F(x). (5) 
In [Su 1975], a partial solution of equations (4) and (5) was also found. 
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The general solution of equations (4) and (5) was found in [Ca 1975] (for 
review of this and the related problems, see also [OP 1983]). This solution 
has the form 

f(x) = a((x;g 2 ,g 3 ) + (3x, (6) 

where ((x) is the Weierstrass zeta-function (see, for instance, [WW 1927]). 

In the present paper, we consider only the three-body problem but in the 
general case when all three particles are different from each other. 

In this case, the ground-state wave function has the form 

ty (x 1 ,x 2 ,x 3 ) = i)i{x 2 -x 3 )ip 2 (x 3 -x 1 )ip 3 (x 1 -x 2 ) (7) 

and satisfies the Schrodinger equation 

-A^o + Uip = E q , (8) 
U = ui(x 2 - x 3 ) + u 2 (x 3 - X]) + u 3 {x x - x 2 ). (9) 

Substituting \l/o from (7) into (8), we obtain 

*o _1 A*o = U-E = 3(tf(x 2 -x 3 ) + f 2 2 (x 3 -x 1 ) + f 3 2 (x 1 -x 2 )) 

- (fi(x 2 - x 3 ) + f 2 (x 3 - x{) + f 3 (xi - x 2 )f 

+2 (f[(x 2 - x 3 ) + f 2 {x 3 - an) + f 3 { Xl - x 2 )) ; (10) 

fi = ^-M- 

Hence, for the potential energy U(x 1 ,x 2 ,x 3 ) to have the form of pairwise 
interactions (9), three functions 

f(x) = f 1 (x), g(y) = f2(y), h(z) = f 3 (z) (11) 

must satisfy the functional equation 

(f(x) + g(y) + h{z)f = F(x) + G(y) + H(z), x + y + z = 0. (12) 

From (10) - (12) it results the following expression for the potential en- 
ergies: 

Ul (x) = 3f(x) + 2f'(x)-F(x)+e u 
u 2 (x) = 3g 2 (x) + 2g'(x)-G(x)+e 2: (13) 
u 3 (x) = 3h 2 (x) +2h'(x) -H(x) + e 3 , 
£i + £2 + £3 — Eq 
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2. Let us consider the meromorphic solutions of the equation 

(/(*) + 9iv) + Kz)f = F(x) + G(y) + H(z) (14) 

satisfying the condition x + y + z = 0. 

Let us call the solution of equation (14) nondegenerate if the functions 
f(x), g(x) and h(x) have the pole in a finite domain of complex x-plane. 

The main result of this paper is the following 

Theorem. The general nondegenerate solution of equation (14) in the 
class of meromorphic functions has the form 

(15) 
(16) 
(17) 

(18) 
(19) 
(20) 



/(*) 


= a((x - 


ai;^2,5s) +/3a: + 7i, 


g(x) 


= a((x - 


a 2 ; 92,93) +/3x + 72, 


h(x) 


= a((x - 


0-3; 92, 93) +/3x + 7 3 , 


F(x) 


= a 2 V(x 


- (11,92,93) + 2ia((x 


G(x) 


= a 2 V(x 


- a 2 ;g2,93) + 27"C(^ 


H{x) 


= a 2 V(x 


- 03; 92,93) + 27«C(^ 



± 

3 ; 

3 ; 

7 2 



w/iere 

ai + a 2 + a 3 = 0, 71 + 72 + 73 = 7- (21) 

Proof . The proof of the theorem is divided on several steps. 
Let us begin with 

Lemma 1. The functions (f(x), g{y), h(z)) satisfy equation (14) for the 
corresponding functions (F(x), G(y), H(z)) if and only if the equation 

( f"(x) g"(y) h"(z) \ 
det f'{x) g'(y) h'(z) = (22) 

V 1 1 1 J 

can be solved under condition x + y + z = 0. 

Proof. Let us apply to equation (1) the operator 

d. • — • — 

dy dx' 
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where 



This gives: 

d_ 

dx 

d_d_ 

dy dx 



= A_A 
dx dy ' 



: 2 (f Or) - h'{z)){f{x) + g(y) + h{z)) = F\x) - H'{z\{2?>) 
: 2h"(z)(f(x)+g(y) + h(z)) 

+ 2{f{x)-h\z)){g'{y)-h'{z))=H"{zl (24) 



d-l-.-Oz ■ h"(z)(f(x)-g\y)) 



dy dx 



+ f"(x) (g'(y) - h\z)) + g"{y) (h'(z) - f'(x)) = 0. (25) 



Here we use the fact that d- is a differential operator and that d-h'(z) = 
d-h'\z) = 0. Hence, if functions f(x), g(y), h(z) satisfy equation (1), then 
these functions also satisfy equation (25) which can be obviously rewritten 
in the form (22). 

Conversely, let the functions f(x), g(y), h(z) satisfy (22) and, conse- 
quently, (25). Equation (25) may be rewritten as 

d_ [h"(z)(f(x) + g(y) + h(z)) + (/'(*) - h'(z))(g'(y) - h\z))\ = 0. 

Then there is the function Hi(z) satisfying the following equation: 

h"(z) (/(*) + g(y) + h(z)) + (f (*) - h'(z))(g'(y) - h\z)) = H^z) (26) 

Let us note that equation (26) is equivalent to the equation 

^ [(/'(*) - h'(z))(f(x) + g(y) + h(z))\ = H^z). 

Therefore, there are functions Fi(x) and H 2 (z) such that H 2 (z) = H^z), 
and 

(/'Or) - h'(z))(f(x) + g(y) + h(z)) = F 1 (x) - H 2 (z). (27) 
On the another hand, equation (27) is equivalent to 

^ (/(*) + 9(y) + Kz)f = 2(F 1 (x) - H 2 (z)), 



i.e., there are functions F(x), G(y) and H(z) such that F'(x) = 2Fi(x), 
H'(z) = 2H 2 (z), and 

(/(*) + g(y) + h{z)f = F(x) + G(y) + H{z). 
Thus, Lemma 1 has been proved. 

Lemma 2. Equation (14) is invariant under the following transforma- 
tions: 



/(*) 
F(x) 

g(y) 
G{y) 
Hz) 

H(z) 



f + a ± x + a 2 f(a 3 x + ai), 
F + a 4 x + alF(a 3 x + «i) + 2a 2 cf(a 3 x + «i), 
go + aiV + a 2 g(a 3 y + a 2 ), 
G + a A y + alG(a 3 y + a 2 ) + 2a 2 cg(a 3 y + a 2 ), 
h + a-iz + a 2 h(a 3 z + a 3 ), 
i^o + a 4Z + a\FL{a 3 z + a 3 ) + 2a 2 ch(a 3 z + a 3 ), 

where a k {k = 1, . . . , 4) and c are /ree parameters 

fo + go + h = c, F + G + H = c 2 , «i + a 2 + tt3 = 0. (28) 

This Lemma is proved by a direct calculation. 

Corollary 3. For appropriate values of the parameters a±, a 2 , a 3 , all the 
functions (f(x),g(y),h(z)), (F(x),G(y), H(z)) are regular at x = 0, y = 0, 
z = 0, respectively. 

The proof follows from the fact that the set of poles of the meromorphic 
function of one complex variable is discrete. Thus, we may suppose that all 
the functions are regular at x = 0, y = 0, z = 0. 

Definition 4. Let us call the solution of equation (14) totally degenerate 
if at least one of functions f(x),g(x), and h(x) is linear. 

The next Lemma describes all the degenerate solutions of equation (14). 

Lemma 5. Let (f(x),g(y),h(z)), (F(x),G(y), H(z)) be a totally degen- 
erate solution of equation (14). 

Three cases are possible. 

1. All three functions f(x), g(y) and h(z) are linear. Then 

f(x) = f + fix, F(x) = F + F lX + {h - g 1 )(f 1 - h)x 2 , 
g(y) = go + giy, G(y) = G + G t y + (gi - f\)(gi - h)y 2 , 
h(z) =h + h lZ , H(z) = H + H lZ + (h - gi )(h! - h)z 2 . 
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Here fo, fi,go,gi,ho, hi are free parameters. 
Let fo + go + = c - Then 

Fo+Go+H = c 2 , F l = b+2cf 1 , G 1 = b+2c gi , H 1 = b+2ch 1 , 
and b is a free parameter. 

2. Two of the functions f(x), g(y) and h(z) are linear. For example, 

g(y) = go + giv, Hz) = h o + hz. 

Then f(x) is an arbitrary function and 

g(y) = go+ay, h{z) = h +az, G{y) = G +by, H(z) = H +bz, 

F(x) = [g + h -ax + f{x)f - (G + H - bx). 
Here go, ho, a, b, Go, Ho are free parameters. 

3. Only one of the functions f(x), g(y) and h(z) is linear. For example, 

h(z) = h + h\z. Then 

f(x) = fo + ax + ci exp(Ax), 

F(x) = F + bx + ci exp(Ax) (2c + c\ exp(Ax)), 

g(y) = go + ay + c 2 exp(Ay), 

G(y) = G + by + c 2 exp(Ay) (2c + c 2 exp(Ay)), 
h(z) = h + az, 

H(z) = H + bz + 2ciC2 exp(— A2;). 
Here a, 6, c, Ci, c 2 , A are free parameters, and 

fo + go + h = c, F + G + H = c 2 . 

Proof. 

Case 1. From (22) it follows that f(x), g(y) and h(z) are the arbitrary 
linear functions. The form of the functions F(x), G(y) and H(z) can be 
reconstructed directly from (14) by taking into account the identity 2xy — 
z 2 — x 2 — y 2 . 

Case 2. From (22), we obtain 

f"(x)(g l -h 1 ) = 0. 
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If f"(x) 7^ 0, then g\ = hi, and f(x) is the arbitrary function. The form of 
the functions F(x),G(y), H(z) can be reconstructed immediately. 
Case 3. From (24), we get 

2{f{x)-hi){g\y)-hi)=H"{-x-y). 

If f'(x) and g'(y) are not constants, then according to the classical Cauchy- 
Pexider result [Ca 1821] (see also [Ab 1823]), we obtain 

f(x) -h 1 = c 1 exp(A:r), g'(y) - h = c 2 exp(\x), 

where ci,c 2 and A are free parameters. Therefore, 

f(x) = f + hix + ci exp(Ax), g(y) = g + h x y + c 2 exp(Ay), 

where Ck = c^/A, k — 1,2. The form of the functions F(x),G(y), H(z) can 
be reconstructed easily. Lemma is proved. 

The functions f(x),g(x), h(x) from equation (1) will be regarded as non- 
degenerate solutions of equation (1). 

Lemma 6. For the appropriate values of the parameters fo, go, ho, a\, F , Go 
(see Lemma 2) we have 

/(0) = g(0) = h(0) = 0, h'(0) = 0, F(0) = G(0). (29) 

The proof is easy. 

Lemma 7. For an appropriate choice of the parameters ct\ and a 2 , we 
have the relation f(x) ^ g{x). 

Proof. Suppose on the contrary that 

f{x + ai) - f(ai) = g(x + an) - g{a 2 ) (30) 

for all aii and a 2 in any neighborhood of the point x — 0. Differentiating 
(30), we obtain 

df(x + t*i) = df(x + ai) = ; 

i.e., 

f(x + «i) = f'{ai)x + /(ai). 

contradicting the assumption that the solution is nondegenerate. Lemma is 
proved. 
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Hence, it is sufficient to find all the nondegenerate solutions of equation 
(1) under the following additional conditions: f(x) ^ g{x) and /(0) = g(0) = 
h(0), h'(0) = 0, F(0) = G(0) =0. 

Interchanging x and y in equation (14), we obtain 

(f(y) + g(x) + h(z)) 2 = F(y) + G(x) + H(z). (31) 

Subtracting (31) from (14), we see that 

[(/(:r) - g(x)) - (f(y) - g(y))\ [(f(x) + g(x)) + (f(y) + g(y) + 2/i(z))] 
= (F(s)-G(s))-(F(y)-G(i/)). 

The last equation can be rewritten as 

tp(x + y) = V (x) + V (y) - ^Zf^y ( 32 ) 

where </?(#) = — 2/i(— x), r)(x) = f(x) + g(x), £(x) = f(x) — g(x), j(x) = 
F(x) - G{x), and ip(0) = <//(0) = r?(0) - 7(0) = £(0) = 0, and ip"(x) ^ 0. 

Definition 8. Let us call the solution (</?, 77,^,7) of equation (32) nor- 
malized if the following initial conditions are satisfied: 

£'(0) = 1, t/(0) = 0. 

Lemma 9. ITie map 

(y?, 77, £, 7) -> (</?, 77 + ^2 (7 + h^ 2 )) , (33) 

where h\ and 62 ore parameters and, 62 7^ 0, defines a group action. Each 
orbit of this group contains one and only one solution. 

Proof. The first statement may be checked by a direct computation. To 
prove the second one, let us differentiate equation (32) with respect to y. At 
the point y = we have 

Assuming (p(x) to be regular at x = and (p"(x) 7^ 0, it is easy to check that 
£'(<)) / 0. 
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Applying the transformation (33) with 62 = (£'(0)) _1 > &i = ~~ V(0)/£'(0) 
to the solution (</?, 77, £, 7), we obtain a normalized solution. Lemma is proved. 

In what follows, the solutions are assumed to be normalized, unless the 
contrary is asserted. Let us now consider the functional equation 

(p{x + y) = (p{x) + ip(y) + T(x)r(y)A(x + y), , . 

^(0) = v9'(0) = r(0) = r"(0) = 0, r'(0) = 1. ld4J 

Lemma 10. For any solution (</?, 77, £,7) 0/ equation (32), i/iere an 
unique solution (ip, r, A) 0/ equation (34) snc/i £/ia£ 

r'(x) - 6 3 t(x) 
77(0:) = <p(x)-<p'(x)£(x), (36) 
7 (x) = V(^) 2 , (37) 

where b 3 = £"(0) a /ree parameter. 

Proof. Let (tp,T,A) be a solution of equation (34). Then acting on (34) 
by the operator d- = (d/dx — d/dy), we obtain 

= (f/(x) - <p'(y) + (T'(x)T(y) - T(x)r'(y))A(x + y), 

A(x + y) = - ™, (V (38) 

T'(x)r(y) -t(x) T'(y) 

Hence, we can transform equation (34) to the equation 

tptx + y) = tptx) + tpiy) + r(x)r(y) - - ffij . (39) 

On the other hand, 

r(x)r(y) t(x) r(y) 1 



r(x)r'(y) - T'(x)r(y) t'(x) r'(y) (r(x)/r'(x) - b 3 ) - (r(y)/r / (y) - 63) 

a*) ay) 

where the function £(x) may be expressed in terms of t(x) by formula (35) 
with the free parameter 63. Therefore, 

ip'(x) - ip'(y) 



<p(x +y) = <p(x) + (p(y) + a x )av) 
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Substituting the expressions for i](x) and 7(2;) from (36) and (37), we 
obtain the solution (<£>, 77, £, 7) of equation (32). 

Now let (</?, 77, £,7) be a solution of equation (32). Substituting y = in 
(32), we obtain 

, \ ( \ l( x ) 
<P(X) = V (X) - ^y, 

i.e., 7(2;) = £(x)6(x), where S(x) = r)(x) — (p(x), and our initial conditions 
</(0) = 7/(0) = = (p(0) = 7/(0) are satisfied. 

Hence, 7'(0) = 0, and from the formula for (fi'(x) obtained in the proof of 
Lemma 9, we obtain 

7 ( X ) = -ip'(x)^ 2 (x), T](x) = <p(x) ~ <p'(x)Z(x), 

as asserted in (36) and (37). Let us note that formula (35) may be regarded 
as the differential equation for the function r(x). Solving this equation with 
initial conditions r(0) = 0, r'(0) = 1 we obtain the function r(x). If, more- 
over, we take 63 = £"(0), this function will satisfy the condition r"(0) = 0. 

Substituting now the expressions for £(#), fj{x), j(x) into equation (32), 
we obtain equation (39). 

Let us apply the operator cL to equation (39). We obtain 

d ( ^ = ) = 0. 

\T(x)r'(y) - r'(x)r(y)J 

Thus, we have proved that the functions tp(x) and t(x) determine the func- 
tion A(x) given by expression (38). Lemma is proved. 

So, we have shown how to construct all the solutions of equation (32) 
using the solutions of equation (34). 

Now we describe the general analytical solution of equation (14). 

Lemma 11. Let (tp,r,A) be a solution of equation (34). (Let us recall 
that <p(0) = <//(0) = r(0) = r"(0) = and r'(0) = 1.) Then the function 
u(x) = (fi'(x) is a solution of the equation 

(u'f = c 3 m 3 + 4c 2 m 2 + 2 Cl u + c 2 , «(0) = 0, u'(0) = c , (40) 

and if c = 0, then c\ 7^ 0. 
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The functions t(x) and A(x) satisfy the following equations: 

t'(x) 1 u'(x) + c 

t(x) 2 u(x) 

A'(x) 1 u'(x) - c 

A(x) ' 2 u(x) 



(41) 
(42) 



If c = 0, then u(x) — \c\ r(x) 2 , and A(x) = \c\ r(x). 

Proof. Let us consider the first three derivatives with respect to y of 
equation (34) 

p'(x + y) = (p'(y)+T(x)[T'(y)A(x + y) + T(y)A'(x + y)}, 
<p"(x + y) = <p"(y) + t(x) [r"(y)A(x + y) + 2r'(y) A'(x + y) 

+r(y)A"(x + y)}, 
tp'"(x + y) = V "\y) + r(x)[T"'(y)A(x + y) + 3r"(y)A'(x + y) 
+ 3T'(y)A"(x + y)+T(y)A"'(x + y)}. 

Taking y = and making use of the initial conditions for <p(x) and t(x), 
we get 

(p'(x) = t(x)A(x), (43) 
<p»( x ) = v "(0) + 2r(x)A'(x), (44) 
<p m (x) = cp"\0) + t(x)It"\0)A(x) + 3A"(x)]. (45) 

Let if k = y? (fe) (°) and r 3 = r w (0). From (43) and (44) we obtain 

cp'(x) A(x) ' 1 ; 

From (45) and (43) it follows that 

^"(x)-ip 3 r 3 A(x) + 3A"(x) 



ifi'(x) A(x) 
Making use of the identity 

A" (A'\ (A'V 
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(47) 



for the quantity <p'(x) —u(x), we obtain the following equation (see equations 
(46) and (47)): 

u" - y?3 / 1 u' - ip 2 V 3 fu' - ip 2 \ 2 

This equation may be rewritten as follows: 

4(«" - y? 3 )w = 4 r 3 w 2 + 6 [uu" - «'(u' - <p 2 )] + 3 (u' - y? 2 ) 2 , 
2 W - 3(w') 2 + 4r 3 w 2 + 4y9 3 w + 3y? 2 = 0. 



(48) 



Let 

T3 = C 2 , ^3 = Ci, </? 2 = Co- 

Equation (48) admits the integrating factor u~ A u' and may be reduced to 
the following equation 

{u-\u') 2 )' = Ac 2 (u~ v )' + 2 Cl (u- 2 y + cg(«- 3 )'. (49) 

Integrating (49) and multiplying the result by u 3 , we obtain equation (40), 
where c 3 is the integration constant. Equation (42) follows from (46). Then 
from equation (43) we obtain: 

u'(x) = t'(x)A(x) +t(x)A'(x). 

From (44) it follows that 

t(x) A\x) = ^(u'(x)-c ). 

Making use of this fact, we obtain 

t'(x)A(x) = ^(u'(x) + co). 

Dividing this equation by equation (43), we come to equation (41). Note 
that if c = 0, equations (41), (42), and conditions r(0) = 0, r'(0) = 1 imply 

U(x) = y T{x) 2 , A(x) = y T(x). 

In particular, it follows that C\ 7^ if c = 0. Lemma is proved. 
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Consider the Weierstrass function p(x) with parameters g 2 and g 3 . We 
have 

p\xf = Ap(xf -g 2 p{x) -g 3 . 

Lemma 12. The general solution of equation (40) may be written in one 
of the following equivalent forms: 



u(x) 



u[x) 



c 3 



(p(x + a) ■ p(a)), 



cMx) + ^*P(x) 2 I r (1 r'(,r). 



where 



2 

1 1 



2 p(x) - c 2 /3 ' 
i/ere p(x) is the Weierstrass function with parameters 



S 2 = 3(^) 



C\C 3 

2 ' 
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_ pc 2 y + cic 2 c 3 _ ^c c 3 y 
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and 



p(a) = ^c 2 , p / (a) = I C oC 3 . 



Proof. Formula (50) gives: 

16 



(u'(x)Y 



r 2 
c 3 



4p(x + a) - g 2 p(x + a) - g 3 



On the other hand, 
4 



{n\x)f = c 3 
+ 4c 2 

Hence, 



Lc 3 



(p(x + a) - p(a)) 



4 

c 3 



(p(x + a) - p(a) 



+ 2ci 



(50) 
(51) 

(52) 
(53) 



c 3 



(p(x + a) - p(a)) 



+ 4 



16 [Ap(x + a) 3 - g 2 p(x + a) - g 3 ] 

= 4 3 [p(x + a) - p(a)f + 4 3 c 2 [p(x + a) - p(a)] 2 

+8cic 3 [p(x + a) - p(a)\ + c\c\. 
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Let us compare the coefficients of terms of the same degree in p(x + a). 
This shows that formula (50) with parameters #2 and g% follows from (53). 
To deduce (51) from (50) one makes use of the addition theorem for the 
p-function (see, e.g., [WW 1927]). 

p(x + a) - pia) = - Mx) + 2 P (a)) + \ ( gig"' )' • 
Therefore, 

\2 



(p(x + a) - p(a)) (p(x) - p(a)) 
= ~ (p( x ) + 2p(a)) (p(x) 2 - 2p(x)p(a) + p(a 

+ - (Ap(xf - g 2 p{x) -g 3 - 2p\x)p\a) + p\af) 



= 3p(x)p(ai) -2p(a) - — p(x) - - g 3 - - p'(x) p'(a) + 



Hence, 



p(x + a) - = \ f ^ p 7 (a) 

2 (p(a;) - p(a)) 2 

3p(a) 2 -^ 2 /4 , 1 



+ 75 y - P (aV + i (p,(a)p(g) - p(a))2 - (54) 

This gives: 

p'(«) = ^ c oC 3 , 3p(a) 2 - -02 = g c i c 3- 

Formula (51) follows from equation (54) by dividing by C3/4. Lemma is 
proved. 

Corollary 13. T/ie general solution of equation (40) /10s £/ie /orm 
. , /cosh2 v /c 2 "a; - l\ sinh2 v /c 2 "a; 

"• w = c '( (w ) +c "^eT (55) 

as c 3 — > 0. 
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Proof. Let 

u*(x) = lim u(x), ip*(x) = lim ip(x), p*(x) = lim p(x). 

C3— >1 C3— >0 C3— >0 

By Lemma 12, the function p*(x) satisfies the equation 

(p'Mf = 4pM 3 -z(^) 2 p*(x) + (^ 3 

C2\ 2 ( , s 2 



4 (p*{x) - (p*(x) + -c 2 j. 



Therefore, 



-p'*(x) 



4 \(p*(x) - c 2 /3) 2 
+ 2c 2 /3 



2^*(x) +4c 2 ^*(a;). (56) 



(p*(x) -c 2 /3) 2 
Differentiating (56) with respect to x, we obtain 

^(x) = 4c2^(x) + 1, V*(0) = 0, Vi(0) = 0. 

Therefore, 

cosh2 v /c 2 x — 1 

In view of (51), 

u*(x) = ci ip*(x) + c tp*(x). 

Corollary 13 is proved. 

Note that according to Lemma 11, if the functions (</?, r, A) satisfy equa- 
tion (34), then the function t(x) is detrmined uniquely by the equation 

r'(x) 1 u'(x) + c 
t(x) 2 u(x) 

subject to the initial conditions r(0) = 0,r'(0) = 1, and the function A(x) is 
determined by equation (43): 

m = 

r(x) 
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Hence, we may regard the functions <p(x) as solutions of equation (34). 
Theorem 14. The general solution of equation (34) 

tp(x + y) = ip(x) + tp(y) + T(x)r(y)A(x + y) 

is given by the function 

V {x) = - (((a) - <Z(x + a)- p(a)x), <p(0) = V '(0) = 0, (57) 

where ((x) and p(x) are the Weierstrass (- function and p-function with pa- 
rameters #2 and g 3 (see Lemma 12). 

Proof. According to Lemmas (11) and (12), it is sufficient to prove that 
any function (p(x) given by formula (57) is a solution of equation (14). It is 
convenient to consider two different cases. 



Case 1. c 3 = 0. 



(pjx) = lim (p(x). 

C3^0 



In this case, (p*(x) = J °° u*(x)dx and hence, using Corollary 13, we 
obtain 

. , sinh2 v /c 2 "a; - 2 v /c 2 "x cosh 2 v /c 2 ~x - 1 

v - {x)=Cl — <w — +c ° (W ■ (58) 

Using the elementary identity 

(59) 

we obtain 

sinh 2y / c 2 "(a; + y) = sinh 2 A /c 2 ~x + sinh l^fc^y 

+ 4 sinh a/c-^e sinh y^y sinh (x + y) , 

cosh 2a/c2"(x + y) = cosh 2^/c^x + cosh 2^/c^y 

+ 4 sinh a/c2~£ sinh cosh a/c^o; + y) . 

Hence, 

V?*(x + y) = (p*(x) + tp*(y) + n(x)n(y)A*(x + y), 

where 

sinh^/c^r ci sinh 

r*(x) = -= — , A*(x) = — -= h c cosh v^x. (60) 

V°2 2 ^ 
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Case 2. C3 7^ 0. Then without any restriction we may take C3 = 2. Accord- 
ing to the Frobenius-Stickelberger formula [FS 1880], the functions 
f(x),g(y),h(z) constitute a solution of equation (1): 

f(x) = C( a i — a/2 — x) — p(a)x — C( a i ~ a /2), (61) 
g(y) = ((-a 1 -a/2-y)-p(a)y + ((a 1 + a/2), (62) 
h(z) = ((a - z) - p{a)z - ((a). (63) 

Using the reduction of (1) to equation (14) described above, we obtain 

ip(x) = - 2h(-x) = 2 (((a) - ((x + a) - p{a)x) 

which gives the solution of equation (34). Theorem is proved. 

Corollary 15. The general normalized solution of equation (32) is given 
by the formulas 

4 

<P(x) = —(C(a)-C(x + a)-'y(a)x), 
c 3 

(64) 

*( X ) = ^ 

' c -2b 3 u(x)+u'(xy 

where 

4 

u(x) = ip'(x) = — (p(x + a) - p(a)) 
c 3 

and 63 is a free parameter, 

r](x) = <p(x) - <p'(x)£(x), 7(x) = -ip'(x)£(x) 2 . 

The proof follows from Theorem 14, formula (61), and from Lemma 10. 
Let us recall that in the proof of Lemma 10 we gave an explicit construction 
of the solution to equation (12) using the solution of equation (14). 

Thus, it is already proved that if (f(x),g(y),h(z)) is the nondegenerate 
solution of equation (1) satisfying the additional conditions 

f(x) = g(x), /(0) = g(0) = h(0) = h'(0), (65) 

then it is necessary to have 

2 

h(x) = — (C(a-x) -~f(a)x-((a)), (66) 
c 3 
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where c^,a and the parameters g 2 , <?3 of the p-Weierstrass function satisfy 
the condition of Lemma 12. Moreover, if C3 7^ 0, then for the functions 

f(x) = —(c(ai-^-x^-p(a)x-c(ai-^yj, (67) 

g(x) = —(c(-a!-^-x^-p(a)x + c(ai + ^y), (68) 

where a is the free parameter, the function h(x) of form (66) gives the solution 
of equation (1). Hence, there are two unsolved problems. 

1. Are the functions f(x) and g(x) for C3 7^ the only functions that give 

the solution of equation (1) for a fixed function h(x)7 

2. How can we find sufficient conditions for C3 = on the parameters of the 

function h*(x) = lim C3 ^o h(x) such that there exist functions f(x) and 
g(x) for which (f(x),g(x),h*(x)) is the solution of equation (1) and 
how can we find all such functions (f(x), g(x))7 

Let us note that in the case C3 = the main problem is that we cannot 
pass to the limit as c 3 — > in formulas (67), (68) (in contrast to (66)). 

To solve these two problems we shall first consider the reduction of equa- 
tion (1) to equation (12) and shall use the general analytic solution of equa- 
tion (12) (see Lemma 9 and Corollary 15). 

Let us begin with the case C3 7^ 0. 

Lemma 16. Let the functions (fi(x),gi(x),hi(x)) satisfy equation (1) 
and the initial conditions under consideration. If hi(x) = H{x) is the func- 
tion from equation (66), then 

A 0*0 = Sif(x) + s 2 g(x), (69) 
9l (x) = t 1 f(x)+t 2 g(x), (70) 

where f(x) and g(x) are given by equations (67) and (68), and si + s 2 = 

i,h + t 2 = 1. 

Proof. For the functions given by equations (47) and (48), we have 
£0*0 = f(x)-g(x) 

C («i - 1 - + C («i + 1 + ^ 



c 3 

C («l - 77 ) - C [Oil + ^ 



(71) 
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Then 



c 3 
2 

c 3 



p cti - — - x 1 -p(a 1 + — + x 



a 



a 

V \ax---x 



a 



P [oii + - + x 



We see that if the parameters a and a\ are sufficiently close to the point 
x = 0, then £'(0) 7^ 0, and the value {"(0) gives the value of free parameter 
63 required to construct the general normalized solution of equation (32). 
Therefore, in this case the general solution of the equation has the form 

tp(x) = -2h(-x), rj(x) b 2 £(x), 

where h(x) is the function (66), £(x) = f(x) — g(x) and rj(x) = f(x) + g(x) 
for the functions (47) and (48). 
Now if we introduce 

fi(x) + gi(x) = r)(x) + b^(x), fi(x) - g x (x) = b 2 £(x), 

then we get 

A 0*0 = ^ r]{x) + bl I b2 £(x) = sxf(x) + siff(g), 



0i (a) 



i ?7(x) + &1 &2 = *i/(a;) + * 2 ^(a;), 



where 



1 , 61 + 



S2= 2 



&1 + 6 2 1 &1 ~ &2 

2 ' 1_ 2 + 2 ; 



t 2 = - 



61 - 6 2 



Lemma is proved. 

Now it remains to find the values of parameters s± and ti for which the set 
of functions (fi(x), g(x), h(x)) from Lemma 16 gives the solution of equation 
(!)• 

Let us introduce the notation 



f"(x) g"(y) h"{z) 
det(/,0,*O = det| f{x) g'{y) h'(z) 

1 1 1 
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and use the following formula (see [WW 1927], p. 458) 

1 o-(x + y + g) ajx - y) <r(y - z) a(z - x) 
- det(p(x), p(y), p(z)) = . 

2 (j d (a;) cr 3 (y) a A {z) 

If the conditions of Lemma 16 are satisfied, we have 

det(«i/(x) + s 2 g(x), hf(x) + t 2 g(x), h(z)) 
= si*i det(/(x), f(y), h(z)) + s 2 t 2 det(g(x), g(y), h{z)) (72) 

On the other hand, 

^ det( f(x)J(y),h(z)) 
c\ , , { { a \ ( a 



= -g- det [p [on - - - xj , p \ai - - - yj , p(a - 2) 

c| a(2ai)a(y — x)a(z — y + a-y — 3a/2) cr(:r — z + 3a/2 — «i) 
4 o" 3 (ax — a/2 — x) <x 3 (ai — a/2 — y) <x 3 (a — 2) 

(73) 



^ det(^),^(y),/i(z)) 

= ^ det(p ^-ai - ^ - x^j , p ^-ai - ^ - , p(a x - 2) 

c| o-(2ai) <r(y - x) cr(y - z + ag + 3a/2) a(x - z + a x + 3a/2) 
4 a 3 (ai + a/2 + x) a 3 (ai + a/2 + y) <r 3 (a - z) 

(74) 

The comparison of expressions (73) and (74) shows that if a.\ = Uk is 
the one of the three halfperiods of the Weierstrass-function p(x), then det (•) 
given by formula (72) is equal to zero identically for all values of Si and t±. 
However, if 011 ^ Uk, k — 1, 2, 3, then this determinant is equal to zero if and 
only if siii = 52^2 = 0. So, we have proved our main result. 

Now let us consider the case c 2 — > 0. 

In this case, the general normalized solution is given by function (58). 
Let us denote 

<p**(x) = lim (p*(x), r**(x) = lim r*(x), A* m (x) = lim A*(x). 

C2^0 
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From (58) we obtain 
According to formulas (40), we have 



x 3 x 2 



= Cl 3[ + c °y- ( 75 ) 



t**(x)=x, A**(x) = Ci ^ + c . 



Further, 



U*) = 7^ ( 76 ) 
V~(x) = ^ + ^-[^ + ^1^-, (77) 

7«(x) = -( C iy + ^)(T^)2- ( 78 ) 

Hence, in this case the general solution of equation (32) is given by the 
functions 

(/2**(x), + ^2^*** 

Corollary 17. T/ie general solution of equation (14) m the class of entire 
functions has the form 

f(x) = «i e Ax + /3x + 71, F(x) = ( ai e Xx + 7/v / 3) 2 + 2a 2 a 3 e- A:c , 
^(x) = a 2 e Xx + (3x + 72 , = {a 2 e Xx + 7/v^) 2 + 2a x a 3 e~ Xx , 

h(x) = a 3 e Xx + [3x + 7s , H(x) = (a 3 e Xx + 7/v^) 2 + 2a x a 2 e~ Xx , 

w/iere 7 = 71 + 72 + 73- 

For £/ie case m which A — > and /or £/ie corresponding a*,, A; = 1, 2, 3, (3, 
we obtain the solution 

\x + 71, F(x) = 2a 2 (x - ai) 4 + 2aj(x - a^ 2 + 7 2 /3, 
< 2 x + 72, = 2a 2 (x - a 2 ) 4 + 20:7(0; - a 2 ) 2 + 7 2 /3, 

+ 73, G{x) = 2a 2 (x - a 3 ) 4 + 207(0; - a 3 ) 2 + 7 2 A 



(A + A - 2A), a 2 = ^- (A + A - 2&), 

(A + A -2 A), 7 = 7i +72 + 73 - 7^ (A 2 + Pi + A 2 )- 



/(*) = 


ax 2 


o(x) = 


ax 2 


h(x) = 


ax 2 


where 






1 


ai = 




6a 




1 


a 3 = 




6a 
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Appendix 

It is interesting to note that the general solution of the functional equation 
(1) has found the applications in another physical context. In [BB 1994] the 
Lax representation for the system of equations 

Qj = J2( a + + b< lk) Vjk(qj - Qk), j = l,...,n, (79) 

was constructed. This system describes the motion of n particles on the line. 
The particular cases of this system are integrable relativistic (b ^ 0) and 
nonrelativistic (b = 0) Calogero-Moser systems, as well as the Toda systems. 
Within the framework of this paper and following [BB 1994], let us consider 
in more detail the case of three particles. 
For the system of equations 

Qj = J2 V Jk((lj ~Qk), .7 = 1,2, 3, 

k+j 

let us search for the Lax representation L = [L, M] in the form 

L(q) =q d + A, M(q) = (Br) d + C, 

where L, M, A, B, C are (3 x 3)-matrices and A = (Ajj,(qj — qk)), Ajj = 0, 
C = {Cjkiqj - q k )), Cjj = 0, B = (B jk (qj - q k )), r = (1, 1, 1), and q d , (Br) d 
are diagonal matrices whose diagonals contain the coordinates of the vectors 
q = 52, 93) and Br, respectively. 

The Lax representation leads to the equation 

q d + [q d , A'] = [A, (Br) d ] + [A, C] + [q d , C], (80) 

where A' = (A' jk (qj — q k )). Therefore, 

T,Vjk(q J -q k ) = [A,C] JJ , (81) 

k+j 

([q d , C-A'} + [A, (Br) d + C]) jk = 0. (82) 
From (82) we obtain C = A' and hence, by (81), 

V jk (qj - q k ) = A jk A' kj - A jk A' kj = -V kj (q k - qj). (83) 



23 



Further, (82) yields 

3 

^AjkiBji - B kl ) + A! jt A lk - AjiA' lk = 0. 
i=i 

Now let us set 

®jk = (A'jAk - A 3l A' lk )/A jk . 
A direct verification shows that (84) implies the condition 

<5>jk + $fcm + $mj = 0. 

Let us introduce the functions 

h(x) = -A 23 (x) A 32 (-x), b 2 {y) = -A 31 (y) A 13 (-y), 
b 3 (z) = -A 12 (z) A 21 (-z). 

Taking x — q 2 — q 3 , y — q 3 — qi, z — q\ — q 2 and using the condition 

$21 + $13 + $32 = 0, 

we immediately obtain 

b 2 (y)A 23 A' 32 - 6i(x)A3iAi 3 + b 3 (z)A 31 A' 13 
-b 2 (y)A 12 A' 21 + b!(x)A 12 A' 21 - b 3 {z)A 23 A' 32 = 0. 

Similarly, from the condition 

$21 + $23 + $31 = 

we get 

b[(x) = A 13 A' 31 - b 2 (y)A 32 A' 23 + b 3 (z)A 32 A' 23 
-b!(x)A 21 A' 12 + b 2 (y)A 21 A' l2 - b 3 (z)A l3 A' 31 = 0. 

Taking into account the relation 

b[(x) = A 23 A' 32 - A 32 A' 23 , b' 2 (y) = A 31 A' 13 - A 13 A' 31 , 
b' 3 (z) = A 12 A' 21 - A 21 A' 12 , 

and adding equations (87) and (88), we obtain 

b 2 (y)b[(x) -h{x)b' 2 {y) + b 3 (z)b' 2 (y) -b 2 {y)b' 3 {z) +b 1 (x)b' 3 (z) -b 3 (z)b[(x) 
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Equation (90) may be rewritten as 



det 



fb'^x) b' 2 (y) b' 3 (z) 

bi(x) b 2 (y) b 3 (z) | =0. (91) 

V i i 



Recall that in our notation x + y + z = and therefore, by Lemma 1, the 
functions 

f(x) = J b 1 (x)dx, g(y) = J b 2 (y)dy, h(z) = J b 3 (z)dz (92) 

satisfy equation (1) (which is identical to (14)) for the corresponding func- 
tions F(x), G(y), and H(z). 

Thus, we obtain the following result. 

Theorem. The system of equations 

Qj = J2 V Jk((lj ~(lk), j = 1, 2, 3, 

has the Lax representation L = [L, M] of the form indicated above if and only 
if 

Vjk(Qj - 5fc) = - ap' (qj - q k + Xj - A fe ). (93) 

Proof. Suppose the system of equations does have the indicated Lax 
representation. Then, by (83), (91), and (92), 

V 23 (x) = f"(x), V 31 (y) = g"(y), V 12 (z) = h"(z), 

and therefore by the main theorem 

V jk (qj -q k ) = - ap'iqj ~qk~qi), I (j, k). 

In view of condition (21), i.e., a x + a 2 + a 3 = 0, we see that the a ; 's may be 
presented in the form Xj — X k , which means that we have (92). 

The proof of the converse statement is in fact contained in [Ca 1975b]. 
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